Ballistic Hot Electron Transport in Graphene 
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We theoretically study the inelastic scattering rate and the carrier mean free path for energetic 
hot electrons in graphene, including both electron-electron and electron-phonon interactions. Tak- 
ing account of optical phonon emission and electron-electron scattering, we find that the inelastic 
scattering time r ~ 10^^ — 10~^ ps and the mean free path I ^ 10 — 10^ nm for electron densities 
n = 10^^ — lO^'^cm"^. In particular, we find that the mean free path exhibits a finite jump at 
the phonon energy 200 meV due to electron-phonon interaction. Our results are directly applicable 
to device structures where ballistic transport is relevant with inelastic scattering dominating over 
elastic scattering. 



The existence [J of gated two-dimensional (2D) 
graphene layers, where carrier transport controlled by 
an external gate has become possible [2|], provides the 
exciting possibility of novel high-speed electronic device 
structures Q utilizing the high graphene carrier mobil- 
ity 0, 01 • Such fast graphene devices would work in the 
ballistic transport regime, where carrier mobility limited 
by elastic scattering, is essentially irrelevant (i.e. Ic > I, 
with Ic, I, being respectively the elastic and the inelastic 
carrier mean free path), and what matters is the inelastic 
scattering due to electron-electron and electron-phonon 
interactions. Such ballistic devices for ultrafast appli- 
cations can only work if the relevant device dimensions 
are smaller than the inelastic mean free path I, and the 
speed of this device is limited by the inelastic scattering 
time T (< Te, where Tq is the elastic relaxation time). 
In currently available high-mobility (> 20, 000 cm^/Vs) 
graphene samples, Ic{tc) ^ lO^nm(lps), and therefore, 
inelastic scattering will dominate device operations for 
length (time) scales below lO'^nm(lps). 

In this Letter, we calculate the inelastic mean free path 
(l) and the corresponding inelastic scattering rate 
in graphene limited by electron-electron and electron- 
phonon interactions. Our study is motivated by electron 
transport in the hot electron transistor device structure 
which is so designed as to allow electrons to traverse 
the base region ballistically. In such a device scheme, 
highly energetic electrons are injected in the emitter re- 
gion which then travel through the base region ballisti- 
cally before reaching the collector region. The fraction 
of electrons a that reach the collector goes as a ~ e"''/', 
where d is the width of the base region. The mean free 
path is given by I — vt, where v is the Fermi velocity 
of electron and r is the inelastic scattering time, which, 
in general, is a strong function of the injected energy of 
electrons and the electron density in the base region. 

We consider two principal mechanisms contributing to 
inelastic scattering arising from many-body interactions: 
(1) absorption or emission of optical phonons by electrons 
due to electron-phonon (c-ph) interaction; (2) exchange- 
correlation effects induced by electron-electron (e-e) in- 
teraction. For the e-ph interaction, we take into account 
the most dominant phonon mode in graphene - the LO 



phonons at the Brillouin zone center F. This mode shows 
up as the 'G peak' resonance observed in Raman scatter- 
ing experiments with a phonon energy wq ~ 200 meV 
!5] . Acoustic phonons couple very weakly to electrons in 
graphene, and the associated scattering rates are on the 
order of 10^^ even at room temperature Q and can 
be ignored when compared to the scattering mecanism 
discussed above. We also neglect electron-impurity scat- 
tering in our calculation because the mean free paths of 
currently available high mobility graphene samples are 
much longer than the inelastic mean free paths to be cal- 
culated in this paper. 

The inelastic scattering of electrons causes damping 
(i.e., decay) of the quasiparticle state, and the inelastic 
scattering lifetime T{k) is given by the imaginary part of 
the self-energy lmY,{k,E) evaluated on the energy shell 
E = ^k, where £,k = £k — l^ (£fe = fi-vk is the electron 



kinetic energy and 



10 ms is the Fermi velocity 



which is constant in graphene irrespective of electron 
density) is the single-particle energy rendered from the 
chemical potential /i. The inelastic scattering rate (or 
damping rate) 1/r consists of two contributions, given 
by h/2T{k) = ImER_ph(fc,a) + ImSR_,(/c, 6), where the 
first term denotes the contribution to the self-energy due 
to electron-LO phonon interaction 7], 

ImI]|Lph(fc,a) = -^5' E E{"^^(^fc'>^) 

A=±l k' 

S{£,k - 6'A + t^o) + [1 - "f(6'a)]<5(6 - S.k'\ - (^q)} 
1 - Acos(0fc/ - 2(j)k'-k) , . 

2 ' ^ ^ 

and the second term on the right denotes the contribution 
to the self-energy due to e-e Coulomb interaction [H , 



ImE«_,(fc,a) = 

E Y^^^Bi^k'x - 6) + nj.(a'A)]Vk'-k 



(2) 



A=±l h 

Im 



1 + Acos(0fc/ 



e(k'-k,6'A-a + *0+) 

Here, nB{x),np{x) = l/(exp(x/fcBT) =F 1 are the Bose 
and Fermi distribution functions, respectively. For the e- 
ph contribution Eq. g = -{[ihv / a^)./lh/2N McLOq) 
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is the e-ph coupling constant Q, where /? ~ 2 is a di- 
mensionless constant characterizing the rate of change 
of the nearest-neighbor hopping energy with respect to 
bond length, a — 1.42 A is the equilibrium bond length, 
is the number of unit cells, Mc — 2.2 x 10'' mc is the 
mass of a carbon atom (mo is the electron mass). For 
the Coulomb contribution Eq. Vq — 2ne'^ /q is the 
bare Coulomb interaction, e(q,u}) = 1 — Vqli((i,uj) is the 
dielectric function given, within the random phase ar 
proximation, by the electron polarizability n(q, w) 
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FIG. 1: (a) Inelastic scattering rate due to e-ph interaction 
versus the single-particle energy of hot electron E = for 
different electron densities n. n = 1,2,5, 10 x lO'^^cm"^ are 
indicated by the red, green, blue and black lines, respectively, 
(b) Corresponding inelastic mean free path /. 

In the following, we calculate 1/T{k) at zero temper- 
ature, which is a very good approximation for graphene 
since at the usual doping density n = 10^^ — lO'^'^cm"^, 
the corresponding Fermi temperature Tp — 1400 — 4300K 
is much higher than room temperature. First we calcu- 
late the self-energy correction due to e-ph interaction, 
Eq. ([T|). Fig. [TJa) shows the scattering rate l/r ver- 
sus electron energy E — for different values of elec- 
tron density. The gap from to 0.2 eV (which is the 
LO phonon energy ojq) is a characteristic feature of the 
LO phonon absorption process; it results from the Pauli 
blocking by those electrons located within an amount 



of energy uq of the Fermi level, so that decay by elec- 
trons with energy £ [—ujq,uiq\ is forbidden due to the 
restricted phase space. Beyond the gap, l/r behaves 
linearly as k due to the linear dependence on momen- 
tum of the graphene density of states, v{k) — ik/irhv. 
As a result of the gap in l/r, the calculated mean free 
path / = VT [Fig. [TJb)] is infinite within the range of 
the phonon energies, before falling off as ^ 1/E with 
electron energy E. Fig. [T] also shows that the scattering 
rate (mean free path) increases (decreases) with electron 
density, as scattering events become more frequent with 
increasing number of electrons. 




0.2 0.4 0.6 0. 

E(eV) 




'0 0.2 0.4 0.6 0.8 1 

E (eV) 

FIG. 2: Calculated scattering rate (a) and the corresponding 
mean free path (b) of hot electron as a function of energy 
E — for different carrier densities n = 1, 2, 5, 10x10^^ 
cm~^. The inset in (b) shows the calculated damping rate 
r(fc) = ft/2r(fc) scaled by Fermi energy as a function of energy 
divided by Fermi energy for different coupling constant = 
0.8 and 0.4, which correspond to graphene on Si02 and SiC 
substrates, respectively. Note the scaled damping rates are 
independent on the density. 

For e-e Coulomb interaction, we show the calculated 
inelastic scattering rate Eq. ^ and the corresponding 
inelastic mean free path, ^ = wt in Fig. [51 The strength 
of Coulomb interaction is characterized by the dimen- 
sionless coupling parameter Ts = e^/nHv, where k is the 
effective background dielectric constant of the substrate 
[l0| . In the inset of Fig. [D^b), the calculated damping 
rates are shown for different coupling constants = 0.8 
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and Ta = 0.4. Before discussing the scattering time we 
consider the scattering rate l/r due to the e-e interac- 
tion. In conventional parabohc-band semiconductors, an 
electron injected with sufficient kinetic energy can de- 
cay via both plasmon emissions and single-particle intra- 
band excitations (i.e. Landau damping) In doped 

graphene, however, injected electrons cannot decay via 
plasmon emission due to phase space restrictions [l^ . 
Multiparticle excitations, which are excluded in the ap- 
proximations used here, will constitute finite damping of 
the quasiparticles, but the effects of such higher-order 
processes are relatively small in graphene. (Note that 
in undoped graphene even single-particle excitations are 
forbidden, so that the scattering rate within the Born 
approximation is zero due to electron-electron Coulomb 
interaction at T = @.) Since only single-particle ex- 
citations give rise to damping of the quasiparticles, for 
electron energy close to the Fermi level the calculated 
scattering rate due to e-e interaction in graphene, similar 
to the case of 2D parabolic-band semiconductors, is given 
by l/r ~ |efc-^^i=^|^ln|£fc-£'F| Sill- Farther away from 
Ep, however, the dependences of l/r(fc) on k in graphene 
and in parabolic-band semiconductors are qualitatively 
different because both plasmon emissions and interband 
processes are absent in graphene. On the other hand, 
the only independent parameters relevant for the quasi- 
particle scattering rate under the Born approximation 
at T = are the Fermi energy and the dimensionless 
coupling constant r^. Therefore the calculated scatter- 
ing rate due to e-e interaction, which has units of energy, 
must be proportional to Ep, and must be a function only 
of Ek/Ep = k/kp and [see the inset of Fig. [UJb)]. 

In Fig. m^a) the inverse lifetime in units of inverse of 
ps is shown as a function of energy £^ = of the hot 
electrons. Below 200 meV, which corresponds to the op- 
tical phonon energy of graphene, the calculated scatter- 
ing time is only weakly dependent on the carrier density 
of the system. However, for very energetic hot electron 
E > leV the scattering time shows a strong density de- 
pendence and increases by almost a factor of two as the 
density decreases from n = 10^^ cm~^ to n = 10^^ cm~^. 
Our results are consistent with recent experiments of de- 
cay time in ultrafast carrier [l3l|, in which the measured 



decay times are in the 0.07 — 0.12 ps range. The corre- 
sponding inelastic mean free path is shown in Fig. [2Kb) 
as a function of energy. We find that the characteristic 
mean free path of the hot electron with 100 meV energy 
above the Fermi energy is about 100 nm. Again below 
200 meV the calculated mean free path is weakly depen- 
dent on density. 

Combining the results from the e-ph (Fig. 1) and e-e 
interactions (Fig. 2), the total scattering rate becomes 
non-zero below the phonon energy 200 meV, and the 
mean free path / ^ 10 — 10^ nm for injected energy 
E < 200 meV whereas Z ~ 10 nm for E > 200 meV. 

We now propose an interesting device principle based 
on our results in this paper, which are peculiar to the 
many-body effects in graphene. In doped (or gated) 
graphene, the dominant inelastic scattering process of hot 
electrons below 200 meV comes from intraband single- 
particle excitation due to screened electron-electron in- 
teraction; above 200 meV, inelastic scattering due to e-ph 
interaction sets in as electrons are now able to emit LO 
phonons. A lateral hot-electron transistor (LHET) de- 
vice 2M where the electrons travel in the graphene sheet 



of the base region can in principle be fabricated, whose 
operation makes use of the abrupt change in the inelastic 
mean free path due to electron-coupled mode scatterings 
of the injected electrons. Thus, by varying the inelastic 
scattering rate through changing the injection energy, one 
can achieve a significant change in the electron mean free 
path and hence the emitter-collector current. It is also 
possible to use the peculiar scattering rate of undoped 
(or ungated) graphene, which is totally suppressed due 
to phase space restrictions [1]. With the application of 
gate voltage one can easily tune the Fermi level through 
the electron density, so that the damping process via e-e 
or/and e-ph interaction can be activated and deactivated. 

In conclusion, we have calculated the inelastic scat- 
tering rate l/r and the inelastic mean free path I in 
graphene. We find that r ~ 10~^ — 10~^ps and / ~ 
10 — 10^ nm, with a finite step jump at 200 meV at the 
LO phonon energy. Our results have direct relevance to 
ballistic transport in graphene fast device structures. 

This work is supported by US-ONR, NRI-SWAN and 
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